Unconventional scanning tunneling conductance spectra for graphene 
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We compute the tunneling conductance of graphene as measured by a scanning tunneling mi- 
croscope (STM) with a normal/superconducting tip. We demonstrate that for undoped graphene 
with zero Fermi energy, the first derivative of the tunneling conductance with respect to the applied 
voltage is proportional to the density of states of the STM tip. We also show that the shape of 
the STM spectra for graphene doped with impurities depends qualitatively on the position of the 
impurity atom in the graphene matrix and relate this unconventional phenomenon to the pseudopsin 
symmetry of the Dirac quasiparticles in graphene. We suggest experiments to test our theory. 
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I. INTRODUCTION 

The low energy quasiparticles of graphene around K 
and K' Fermi points have Dirac-type properties^. In 
particular, the pseudospin of these quasiparticles around 
K{K') points along (opposite to) their direction of mo- 
tion. The presence of such Dirac-type quasiparticles 
with definite helicity leads to a number of unusual elec- 
tronic properties in graphene^"—. Recently, the influence 
of such Dirac quasiparticles on properties of graphene 
doped with magnetic/non- magnetic impurities have at- 
tracted theoretical and experimental attention^"—. How- 
ever, the recent experimental observation of dependence 
of STM tunneling spectra on the position of the impu- 
rity in the graphene matrix in Ref. lacks a theoretical 
explanation even at a qualitative leveli^. 

Scanning tunneling microscopes (STM) are extremely 
useful probes for studying properties of two or quasi-two 
dimensional materialsii*i^. Studying electronic proper- 
ties of a sample with STM typically involves measure- 
ment of the tunneling conductance G{V) for a given ap- 
plied voltage V . The tunneling conductances measured 
in these experiments have also been studied theoretically 
for conventional metallic systems and are known to ex- 
hibit Fano resonances at zero bias voltage in the presence 
of impuritiesi^ii^. The application and utility of this ex- 
perimental technique, with superconducting STM tips, 
have also been discussed in the literature for conventional 
systems^^. However, tunneling spectroscopy of graphene 
using superconducting STM tips remains to be studied 
both experimentally and theoretically. 

In this work, we compute the STM response of doped 
graphene and demonstrate that the STM spectra has sev- 
eral unconventional features. For undoped graphene with 
Fermi energy Ep = 0, the derivative of the STM tunnel- 
ing conductance (G) with respect to the applied voltage 
[dG/dV) reflects the density of states (DOS) of the STM 
tip [pt), i.e., dG/dV ~ +i-)pt for V > (<)0. By tun- 
ing Ep , one can interpolate between this unconventional 
pt ~ ztdG/dV and the conventional pt ^ G (seen for 
Ep 3> eV) behaviors. Further, we find that for super- 



conducting STM tips with energy gap Aq, G{dG/dV) 
displays a cusp (discontinuity) at eV = —Ep — Aq as 
a signature of the Dirac point which should be exper- 
imentally observable in graphene with small Ep where 
the regime eV > Ep can be easily accessed. For impu- 
rity doped graphene with large Ep, experiments in Ref. 
1^ have seen that the tunneling conductance, as measured 
by a metallic STM tip, depends qualitatively on the posi- 
tion of the impurity in the graphene matrix. For impurity 
atoms atop the hexagon center, the zero-bias tunneling 
conductance shows a peak; for those atop a graphene site, 
it shows a dip. We provide a qualitative theoretical expla- 
nation of this phenomenon and show that this unconven- 
tional behavior is a consequence of conservation/breaking 
of pseudospin symmetry of the Dirac quasiparticles by 
the impurity. We also predict that tuning Ep to zero 
by a gate voltage would not lead to qualitative change 
in shape of the conductance spectra when the impurity 
is atop the hexagon center; for impurity atop a site, the 
tunneling conductance would change from a dip to a peak 
via an antiresonance. 

The organization of the rest of the paper is as follows. 
In SecHIl we present the derivation of the tunneling cur- 
rent. This is followed by Sec. IIIII where we present our 
main results. Finally we conclude in Sec. IIVI 



II. COMPUTATION OF TUNNELING 
CURRENT 

The experimental situation for STM measurement is 
schematically represented in Fig. [TJ The STM tip is 
placed atop the impurity and the tunneling current I 
is measured as a function of applied bias voltage V. The 
possible positions of the impurity is shown in the right 
panel of Fig. [TJ Such a situation can be modeled by the 
well-known Anderson Hamiltoniaiii^. Here we incorpo- 
rate the low-energy Dirac quasiparticles of graphene in 
this Hamiltonian which is given by 

H = Hg + Hd + Ht + HGd + HGt + Hdt (1) 
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The corresponding parameters V°{k), U^.^{k), and W° 
are taken to be independent of valley and spin indices 
of graphene electrons but may depend on their sublat- 
tice index or pseudospin. Note that the tunneling terms 
[Eq. [5] automatically take into account potential scatter- 
ing; such terms are generated once the impurity degree of 
freedom is integrated out unless there is perfect particle- 
hole symmetry {Ep = 0). 

The tunneling current for the present model is given 

by 



FIG. 1: (Color online) Schematic experimental setup with the 
right panel showing two possible positions (atop hexagon cen- 
ter and atop a B site) of the impurity. The numbers denote 
nearest neighbor A and B sublattice sites to the impurity. 
ai(2) = +{—)y/3/2x + 3/2y [lattice spacing set to unity] are 
graphene lattice vectors. The choice of coordinate center (0,0) 
are shown for each case. 
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where N = P^^ii^s is the number operator for the tip 
electrons. These commutators receive contribution from 
Hdt and Hat in Eqs. ([6]) and (Oand can be written as 
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where we define the standard Keldysh Green's function 
Q and G as^ 
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Hat = E (Ua-:,{k)ci,{k)il^+h.c.) (7) These hybrid Green's functions (Eq. HH) obey the usual 

a=Xs;.-^'=^ ' ' ^ T.„,.-„u T. ^(2)< ^(2)> 



Here Hq is the Dirac Hamiltonian for the graphene elec- 
trons which are described by the two component annihi- 
lation operator 



Keldysh relations. For example, Q^J ^ and Q^J ^ can be 
expressed in terms of the time ordered (G^J*), anti-time 

ordered (G^l*), retarded {g!^1^), and advanced (G^l^) 
Keldysh Green's functions aaii 
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belonging to the valley P — K^K' and spin s |, / 
is the identity matrix, r and a denote Pauli matrices in 
valley and pseudospin spaces, vp is the Fermi velocity, 
and = Y.p=K,K'T.s=ni ^'^ denotes the im- 
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Similar relations hold for ci-^j ^{t\ k) and Gu^a ^{t\ k). 
Next, we expand the hybrid Green's functions 



^(i)</ 



purity atom Hamiltonian with an on-site energy and 
U is the strength of on-site Hubbard interaction. Ht is 
the Hamiltonian for the superconducting (Aq 7^ 0) or 
metallic (Aq = 0) tip electrons with on-site energy ct^, 
where v signifies all quantum numbers (except spin) as- 
sociated with the tip electrons. The operators ds and 
t,ys are the annihilation operators for the impurity and 
the tip electrons. The Hamiltonians Hcdi Hat, and Hdt 
describe hopping between the graphene and the impu- 
rity electrons, the graphene and the STM tip electrons, 
and the impurity and the STM tip electrons, respectively. 



G'.'>^{t-k), Gi'-J<{t;k), G'/J^it) and GlV^{t) in pertur- 
bation series"'^^. After some straightforward algebra, one 
obtains, to first order in perturbation theory. 
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(13) 

where all the Green's functions appearing in Eq. 

are at the same time t which we have not written 

out exphcitly for clarity. In Eq. p^ . g^^.j^,^,{t) — 

—i{il^{t)i^:cr'{0)) denotes the Green's function for the 
tip electrons which, in frequency space, is given by 

where f{x) = 1/(1 + exp [x/ZcsT]) denotes the Fermi- 
Dirac distribution function at a temperature T, /i( is 
the chemical potential for the tip electrons, and ks is 
the Boltzman constant. Similar expressions can be ob- 
tained for 5* and 5* using Eq. ([T^ . g<^,{t;k,k') = 

— i(V'CT'(*; ^)'0cr(O; k')) denotes the Green's function of the 
Dirac electrons in the presence of the impurity. The re- 
tarded and advanced components of this Green function 
which we shall need in subsequent analysis can be written 
as 
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where again it is understood that all Green's functions 
are at a given time t and G^^i'^^(i) — —i{dl{t)da-' (0)) 
denotes the retarded(advanced) Green's function of the 
interacting impurity electrons. In frequency space, this 
Green's function is given by 



uj - Ed - Rc[I]d(w)] - (+)ilm[l]d(w)] 

(16) 

where Sd(i^) denotes the self-energy of the impurity in 
the absence of the tip. receives contributions from 



both the on-site Hubbard interaction U of the impurity 
electrons and the coupling of the impurity to the Dirac 
electrons. Note that we have neglected the effect of the 
STM tip while computing t/^j^'' {t; k, k') which is justi- 
fied as long as we restrict ourselves to linear-response the- 

[t; k) denotes the single-particle 



ory. InEq. m.Q^^'^'^^^^ 
Green's function for the graphene electrons in the ab- 
sence of the impurity and the STM tip and is given, in 
frequency space, by 



(cj + Ef)I - hvF{Tz(Txkx + gyfey) 



h^vW^"^ - {+)ir] 



(17) 



Finally, the Green's function G'^^, [t; k) = 
-i{di{t)^p^,{0;k)) used in Eq. (USD is a hybrid Green's 
function whose retarded and advanced components are 
given, within first-order perturbation theory, by 



G'J'^\t;k) = 



Next we follow Ref.[lJto substitute Eqs. (HH) and ([14]) 
in Eq. pHI) and approximate the coupling functions to 
be independent of momentum: C7°(fc) = ?7°, = 
and y°(fc) = Such an approximation is justified as 
long we restrict ourselves to low applied voltages. With 
this approximation, after some algebra involving Eqs. 
([TU| .. (|18I) . one obtains the expression of the current 



Xq I duj[f{uj - eV) - f{uj)]pt{uj ~ eV)^pG{uj) 
, \B{lo)\^ \qiuj)\^-l + 2Re[q{uj)]xiuJ 



x\U"\' 



ImSrf(^) (l+;^2(^))(l+e2) 



(19) 



where Iq = 2e(l + £,'^)/h, pg{^) and pt{e) are the 
graphene and STM tip electron DOS, respectively, ^ = 
|[/0|/|C/0| ^ |Vg|/|V2| is the ratio of coupling of the 
impurity to the electrons in B and A sites of graphene 
with U\ = C/° and Vl ^ and Ed(e) is the impu- 
rity advanced self-energy in the absence of the tip. Here 
B{e) = V^U^l2{f) and q{e) and x(e) are given by 



Xie) 



[W'/U" + V'h{e)]/[V"hie)] 
e~ Ed- ReS]rf(£) 
Imi;rf(e) 



(20) 



where we have neglected the energy dependence of 
the coupling functions assuming small applied voltages. 

In Eq. m. h{e) = (1 +e')E/cTr[Re{gi°)^(e,k)}], 

/2(e) = {l+e) Ek Tr [lm{^^°^'^(e, k)}] , and Tr denotes 
trace over Pauli matrices in pseudospin, valley and spin 
spaces. Substituting Eq. ([12]) in Eq. ^0^, we find^-" 

/i(e) = -4(l + e')(e + SF)ln|l- AV(e + SF)'|/A2 
/2(e) = A{l+^^)TT\e + EF\0{A-e- Ef)/A^ (21) 



4 



where A is the ultraviolet momentum cutoff and 9 is the 
Heaviside step function. Usually, in graphene, A is taken 
to be the energy at which the graphene bands start bend- 
ing rendering the low-energy Dirac theory inapplicable 
and can be estimated to be 1 — 2eV^. 

Equations (|19ll2ip constitute the central results of this 
section and yields the tunneling current through the STM 
tip within linear-response theory. We are going to ana- 
lyze these equations in the subsequent sections. 



III. RESULTS 

In this section, we are going to analyze the tunneling 
conductance {G{V) = dl/dV) as measured by the STM 
tip. First we consider the case of a superconducting tip 
in the absence of any impurity. In this case, the contribu- 
tion to the conductance comes from the first term of Eq. 
(IT^ . For s-wave superconducting tips, one finds that the 
tunneling conductance {G{V) — dl/dV) for Ep > and 
at T = is given by (with r — Ep/Ao, p = —eV/Ao) 



G 

dG 
dV 



Go + / Sgn{z-p + r)Ut{z)dz (22) 



eGo 
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Ut{p)^Ut{p)\r\~20{p-r)Ut{p- 
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where Gq = 87r^e2|C/"p(l-|-g^)po tpoAo/ fe, pa = Po|r-p|, 
Pt{r) = potNt{r), Nt{x) = \x\/^/x^e{\x\ - 1), sgn(x) 
denotes the signum function, po = 6\/3/(27r?i^w|,) (Ref. 
Ili) and Pot is the constant DOS of the metallic tip. For 
graphene with Ep = r = dG/dV - Sgn(y)M(-l/), 
i.e., the tip DOS is given by the derivative of the tun- 
neling conductance. For large Ep away from the Dirac 
point, the first term of G becomes large and reflects the 
tip DOS. In between these extremes, when Ep ^ eV, 
neither G nor dG/dV reflects the DOS. In this region, 
the signature of the Dirac point appears through a cusp 
(discontinuity) in G (dG/dV) at eV = —Ep — Aq arising 
from the contribution of the second (third) term in Eq. 
([^ [Eq. (US])]. These features, shown in Fig. distin- 
guish such graphene STM spectra with their conventional 
counterparts^^. 

Next, we turn to the case of impurity-doped graphene 
and consider a metallic tip with constant DOS. The con- 
tribution to the tunneling conductance from the impurity 
(after subtracting the graphene background) at T = 
[Eq. (in])] is 



Gil 



Gr 



\B{V)\^ \q{V)\^-l + 2Re[q{V)]x{V 



ImSd(F) 



■?4) 



where Gg = 2e^potA/h. Such tunneling conductances 
are known to have peak/antiresonace/dip feature at zero 
bias for |g|^l/~l/<C 1—. In conventional met- 
als or earlier STM studies in graphenei^, U'^ has been 




FIG. 2: (Color online) Plot of the tunneling conductance G 
and its derivative dG/dV as a function of the applied bias 
voltage ey/Ao = —p for r — 0,2,6 (red solid, blue dashed, 
and black dotted lines), respectively. See text for details. 



taken to be a flxed parameter independent of the po- 
sition of the impurity. However, as we show here, the 
situation in graphene necessitates a closer attention. To 
this end, we note that [7° is proportional to the proba- 
bility amplitude of the Dirac quasiparticles in graphene 
to hop to the tip. and its strength can be estimated 
using the well-known Bardeen tunneling formulaF^: 



/ d^r {^Uz)d,^Gir, z) - ^l,{r, z)d,Mz)'^ 

^'^(n), zq), where the last similarity is obtained by a care- 
ful evaluation of the surface integral / d^r over a surface 
between the graphene and the tip parallel to the graphene 
sheet--, (^0, zo) is the coordinate of the tip center, 0,y(z) 
is tip electron wavefunction, and the wave- function of the 
graphene electrons '^'ai'T, z) around K{K') valley, can be 
written, within tight-binding approximation, as^'^ 



G[r,z) 



E 



J[{K(K')+5k}-Rf] 



V{r-Rt) 



+e+(-)*«'=(/j(f-i?f)l/(z) 



(25) 



Here 9k = aictaniky/kx), Sk is the Fermi wave vector 
as measured from the Dirac points with \6k\ ^ |-ftr(if')| 
for all Ep, ifif) are localized pz orbital wave functions, 
iV is a normalization constant, f{z) is a decaying func- 
tion of z with decay length set by work function of 
graphene, and i?^*-^' = ndi -I- md2(a2 ~ y) with inte- 
gers n and m denote coordinates of the graphene lat- 
tice sites (Fig. When the impurity and the STM 
tip is atop the center of the hexagon, pseudospin sym- 
metry necessitates ip{ro — Rf'^) to be identical for all 
neighboring A and B sublattice points 1 — 6 surround- 
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FIG. 3: (Color online) Plot of Gimp as a function of V for 
|W^°/;7"| = 0.05 (right; impurity atop a site) and 2 (left; im- 
purity atop hexagon center) for Ef/A = 0.3, 0.1, andO (black 
solid, blue dashed, and red dotted lines, respectively). Plot 
parameters are 5U ^ V° ^ 0.05A, W° = 0.0005A, and ea = 0. 

ing the impurity (Fig. [ij. Consequently, the sum over 
lattice vectors Rf in Eq. ([25]) reduces to a sum over the 
phase factors exp{i[{K{K') + 6k} ■ Rf]) for these lat- 
tice points. It is easy to check that this sum vanishes 
for both Dirac points (when \6k\ = 0). Thus the only 
contribution to ^g('^Oi-2^o) comes from the second and 
further neighbor sites for which the amplitude of local- 
ized wave functions (p{'ro — Ri ) are small. For finite 
Ep, {5k 7^ 0) there is a finite but small contribution 
(0(|(5fc|/|^|)) to '^airoTZo) from the nearest- neighbor 
sites. Thus ^ai^OiZo) and hence C/° is drastically re- 
duced when the impurity is atop the hexagon center. In 
this case, we expect ^ and hence |g| 3> 1 [Eq. 
(l20t ] leading to a peaked spectra for all Ep. In con- 
trast, for the impurity atom atop a site, there is no such 
symmetry induced cancellation and ■(/'g('\): -zo) receives 
maximal contribution from the nearest graphene site di- 
rectly below the tip. Thus we expect \W^\ (since 
it is easier for the tip electrons to tunnel to delocalized 
graphene band than to a localized impurity level) lead- 
ing to (? ~ h/h ^ -ln|l - A^/{eV + Epyi/TT. For 
large \eV + Ep\ and impurity atop a site, q < 1 leading 
to a dip or an antiresonance in Gimp which is qualita- 
tively distinct from the peaked spectra for impurity atop 
the hexagon center. As Ep ^ 0, q diverges logarithmi- 
cally for small eV. However, it can be shown that in this 
regime x shows a stronger linear divergence for eV ^ 
which suppresses Gimp. At eV = e^, the divergence of x 
also becomes logarithmic and we expect a peak of Gimp. 
Note that these effects are independent of and hence 
of the precise nature of the impurity. Such an impurity 
position-dependent peak/dip structure of Gimp has been 
observed for magnetic impurities in Ref. |^ for Ep 3> eV . 

To demonstrate this feature, we restrict ourselves to 
impurities with small Hubbard U and compute the self- 



energy of the impurity electrons within a mean-field the- 
ory where Un„ng = U{n(j)ng leading to spin-dependent 
on-site impurity energy e„ — ed + U{ngy . Using Eqs. 
(H)) and (IS|), one then obtains the mean-field advanced 
impurity Green's function C/™p(cj) — {uj — e„ — T,d{io))^^ 
where the impurity self-energy is given by Sd(w) = 
|T^°P(/i + i/2) and mean- field self-consistency condition 
demands Ua — J doj / 'KlniQ"'^^ {uj) . Following Ref. i?!, we 
solve these equations to get xi^)^ and ImEd(e) which can 
be substituted in Eq. ([M)) to obtain Gimp- We note, from 
Eqs. ([20]) and (p4| . that Gimp/Gg depends on the ratios 
Ep/A, V"/A, and W°/U° which can not be quantita- 
tively determined from the Dirac- Anderson model. We 
therefore treat them as parameters of the theorjiiii and 
compute Gimp for their representative values as shown 
in Fig. [21 In accordance with earlier discussions, we 
find that for large Ep/A = 0.3, Gimp has qualitatively 
different features; for the impurity at the center of the 
hexagon, it shows a peak (left panel) while for that atop 
a site (right panel), it shows a dip. The change in Gimp 
from a dip to a peak via an antiresonance as a function of 
Ep/A when the impurity is atop a site can be seen from 
right panel of Fig. [31 In contrast, the left panel always 
shows peak spectra. 



IV. CONCLUSION 

In conclusion, we have shown that the tunneling con- 
ductance spectra of both doped and undoped graphene 
have unconventional features not discussed in earlier 
studiesiS. In particular, the STM spectra of doped 
graphene depend qualitatively on the position of the im- 
purity in the graphene matrix. This feature is demon- 
strated to be a direct consequence of pseudopsin symme- 
try and Dirac nature of graphene quasiparticles. 

Further experimental verification of our work would in- 
volve measuring tunneling conductance of doped and un- 
doped graphene by varying Ep. For undoped graphene 
with Ep — 0, we propose to measure the tunneling con- 
ductance spectra using a superconducting tip and ver- 
ify that dG/dV ~ ptSgn{V). For smaU Ep > 0, where 
it is possible to access the regime eV > Ep in experi- 
ments, we predict a cusp (discontinuity) in G{dG/dV) 
at eV — —Ep — Aq as a signature of the Dirac point. 
The variation in the shape of the spectra of impurity- 
doped graphene with impurity atop a site may also be 
experimentally studied. 

We also note that the theory of tunneling conductance 
derived here should also be applicable to the impurity- 
doped Dirac electrons on the surface of strong topological 
insulators with a single Dirac cone^^. In this case, we 
expect to find that the STM spectra should change from 
a dip to a peak through an anti-resonance as the Fermi 
energy is tuned toward the Dirac point. This behavior is 
qualitatively similar to that shown in the right panel of 
Fig. [31 However, such a controlled tuning of Fermi energy 
of topological insulators seems to be experimentally more 



difficult than graphene. 

Note added. Recently, we came to know of Ref. [22 
with related results. 
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